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Abstract 

We analyse the situation where the primordial curvature perturbations are pro- 
duced by the joint effects of an inflaton field and two curvaton fields. We present 
general equations which allow the reader to obtain /jvx f° r several different scenarios 
which differ in the order in which fields decay into radiation after inflation. In order to 
investigate the physics of these equations we analyse some simplified situations where 
the fields are harmonic and both curvatons are frozen at the same expectation value 
during inflation. We find quite complex behaviour - for a given situation where the 
inflaton contributes a fixed amount to the total curvature perturbation there are situ- 
ations where f^L is maximised if both curvatons share equally in contributing the rest 
and situations where fx l is maximised if only one of the two curvatons contributes the 
rest. We are unable therefore to make any completely general extrapolations about the 
expected non-gaussianity from N curvatons. We find that as the curvaton contribu- 
tion to the overall perturbation is gradually increased, Jnl rises to a maximum before 
falling again and that a given /jvl can correspond to many different parameter sets for 
the two curvatons. 

1 Introduction 

Inflation is the most successful theory we have to describe the early universe. According to 
inflation [lj, at some early time the universe underwent an accelerated expansion, and this 
solves the horizon and flatness problem while also predicting an almost scale-invariant power 
spectrum in the cosmic microwave background, consistent with current observations. The 
most basic models of inflation (see, for example, j2j) involve a single, massive scalar field 
called the inflaton, which slowly rolls down its potential to source the necessary expansion. 
Once it reaches the minimum of the potential it oscillates around that point, reheating 
the Universe by decaying into standard model particles. The quantum fluctuations of the 
inflaton get stretched out with the expansion and become classical perturbations, causing 
inhomogeneities which upon horizon re-entry evolve to become the large-scale structure that 
we observe today [3]. 



1 



Since many theories that go beyond the standard model predict scalar fields, it is not 
unreasonable to assume that in early times, these fields had expectation values away from 
the minimum of their potential. Such fields, known as curvaton fields, can be present during 
inflation but only begin to move when the Hubble parameter drops down to their mass (as 
in [3] and [3]). They then proceed to oscillate around the minimum of their potential and 
eventually decay when the Hubble parameter is equal to their decay rate. The presence of 
these particles imprints a record of the curvaton perturbations as well as the inflaton pertur- 
bations upon the cosmic distribution of matter and therefore provides a different mechanism 
for the creation of the primordial curvature perturbation ( and can alter the predictions of 
basic single-field inflation [Tj, [8], [H], [ID]- There are a plethora of curvaton models one 
can imagine existing in particle physics and string theory and they typically include a large 
number of parameters. To constrain them, one needs to compare their predictions to current 
astrophysical observations, with the most popular discriminant being the non-gaussianity 
parameter, f NL ([II], [12],[T3], [TJ]). 

The idea behind the curvaton scenario is to boost small perturbations and make them 
larger by virtue of the energy density in the curvaton field red-shifting more slowly than 
the energy density from the inflaton field. This occurs once the inflaton has decayed into 
radiation but the curvaton is still oscillating (and thus red-shifting like matter). There are 
therefore multiple scenarios depending upon the order in which the inflaton decays, the 
curvaton starts to move and the curvaton decays. In this work we will attempt to study 
a number of these different situations, although we will not consider the case where the 
curvaton starts to roll during inflation, as this moves into the territory of multifield inflation. 

In inflation models, the non-gaussianity parameter f^L is related to the curvature per- 
turbation after the last decay. This parameter is a measure of deviations from a purely 
Gaussian distribution of the primordial curvature perturbation and the current f^L limits 
are given by WMAP 7 to be —10 < f^L < 74 [IS]. We analytically compute the expression 
for C at the last decay and /jvx i n models with one inflaton, <ft, and two curvatons, a and \i 
for various sequences of particle decays. In each case, the inflaton is the one responsible for 
producing the necessary e-folds and both the curvatons are frozen during the inflaton's slow- 
roll. Once inflation has finished, begins to oscillate around the minimum of its potential 
and then decays. Independently, some time after the end of inflation, when H = m curvaton , 
the curvatons begin to oscillate and eventually decay. All three particles contribute to the 
curvature perturbation. We will evaluate the primordial curvature perturbation at the time 
of the last particle's decay, whether it is the inflaton or one of the curvatons. 

2 Perturbations and Non-Gaussianity with multiple fields 

In this section we will present our notation and derive some general expressions for the 
curvature perturbation £ and the non-gaussianity parameter f^L which we will be able to 
use in the following sections. 

2.1 Curvature perturbation from the curvatons 

The simplest potential used in inflationary models is the quadratic, V = |m|</> 2 which we 
will use later on for our examples but in our equations we are allowing the curvaton potential 
to deviate from this case away from its minimum. 
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For each of the two curvatons, the oscillation amplitude will then be a function of the 
initial field expectation value at horizon exit, i.e. a = c(cr*) and x — d{x*) respectively 
where the index * means the expectation value at horizon exit. This general case simplifies 
to the quadratic potential when c = a* =>■ d = 1 and d — x* =^ d! — 1- The curvature 
perturbation up to second order for a is derived in [TH] and it is 



c - 26a (1) 

C<7( 2 ) = ~2 ( X ~ ~d*) ^ 2(l) ^ 

the expressions for the other curvaton, x wm take a similar form. For each case that we 
study we will use C<^>, Co- and ( x to define the isocurvature perturbation, Si oc Q — i = a, x, 
to describe the non-adiabatic part of the curvature perturbation of the curvatons. Next, 
we define the general curvature perturbation in terms of C^, S a and S x to look at the non- 
gaussianity parameter /jvl- 



2.2 General expression for f^i 

In this section, following a similar formalism as in [IB], we derive the general expression for 
the Jnl parameter in terms of the coefficients in the curvature perturbation produced in a 
model with one inflaton and two curvatons, taking into account contributions from all three 
particles. 

Up to second order, the general form of the primordial curvature perturbation will be: 
C = A C/i(i) + BS^i) + TS X{1) 

111 

+ 2^^(1)^(1) + 2^ > ^W^ X W + 2^ ct(1) ^ x(1) 

evaluated at the time of the last decay (the point at which all three fields have decayed into 
radiation). Since /jvl is related to the primordial bi-spectrum, we compute the three point 
correlator of the curvature perturbation: 

(C(k 1 )C(k 2 )C(k 3 )) = {[A 2 AP c ^(k 1 )P Ct M+B 3 EPs 9l (k 1 )P s ^) 

+T 2 ZP Sxl (k 1 )P Sxl (k 2 ) + ^ABHPcMPs^fa) 

+ 1 -AT&P c Jk 1 )P Sxl (k 2 ) + hsTIPs^kjPsJk,)] 
+permutations}(2vr) 3 5 {3) (ki + k 2 + k 3 ) (4) 

where we have assumed that the adiabatic perturbation is uncorrelated with the isocurvature 
perturbations, so (C^(ki)S ff(1) (k 2 )) = (C^ (1) (ki)S X(1) (k 2 )) = {S a(1) (ki)S X(1) (k 2 )) = 0. By 
looking at equation p|, we can now see that the power spectrum for £ at leading order is 
just given by 

P ( = a 2 p u + B 2 P 5ct + T 2 P 5x (5) 
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Where the power spectrum for each component is |lb].|17|. |lUj : 



here e* is the inflaton slow-roll parameter and Mp\ is the reduced Planck mass, Mp\ = 
(87TG)- 1 / 2 . 

Combining equations ^ , (|6| , ^ and ^ we can express all three individual components 
of the power spectrum P^,, Ps a and Ps x in terms of the total P^, and thus substitute them in 
the bi-spectrum equation Q. The non-linearity parameter is defined in the following way: 

(C(k 1 )C(k 2 )C(k 3 )) = B(k u k 2 , fc 3 )(27r) 3 5( 3 )(k 1 + k 2 + k3) 

6 

B(k 1 ,k 2 ,k 3 ) = -f NL [P Ci (ki)P Ci (k2) + permutations]. 



And so we are able to obtain the following equation for the /tvl parameter: 

5 1 2 4 4 

fNL = 6 [A 2 c 2 d 2 + m 2 M 2 Pl e*d 2 d? + 8T 2 M Pl e*c 2 d> 2 ] 2 x ^ A Ac d 

+64B 2 EM^e 2 c'V + MT 2 ZM Pl e 2 ,c 4 d' 4 + 4ABRM Pl e*c 2 c' 2 d 4 
+4ATQM Pl e*c 4 d 2 d /2 + 32BTlM Pl elc 2 'c' 2 l d 2 'd' 2 } (9) 

This complicated looking expression is the general form of the non-linearity parameter 
when one inflaton and two curvatons are all contributing to the total curvature perturbation. 
The particular values of the parameters A, B, T etc. depend upon the particular cosmic 
history. In the following sections we analytically compute the coefficients A, B, T, ... for 
several different cases. However before we do that we will check for consistency with previous 
results. 



2.2.1 First Consistency Check - f^L for curvatons without perturbations from 
inflation. 

Before looking at the relatively complicated situation of mixed perturbations from an inflaton 
and two curvatons, we would like to check that our general equation ^ is able to re-create 
the results of others in particular limits. The first such consistency check, i.e. the case 
where the perturbations come from two curvatons with no contribution from inflation, can 
be done immediately. The other two consistency checks we will make (expression for a single 
curvaton and mixed perturbations from one inflaton and one curvaton) will have to wait 
until the parameters A,B,r, A... etc have been defined in the next section. 

If in our analysis above, we assume that the inflaton's contribution to the primordial 
curvature perturbation can be ignored, we have the case where the perturbations are coming 
from two curvatons alone. In this case, A = A = H = @ = and it is straightforward to see 
that the expression for f NL becomes 

5 B 2 E + ^A 2 Bn + A 4 r 2 Z 
Tnl - g (B2 + A 2 r 2 ) 2 ^ ' 
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where A = cr*/x*- We have also assumed linear evolution of the two curvaton fields since 
horizon exit and so c = a* and d — %*■ This is in agreement with the expression derived in 
[16] for this particular case. 



3 Mixed perturbations from an inflaton and two cur- 
vat ons 

Now we move on to the main part of the analysis in this paper, namely the situation where 
the efolds of inflation required to solve the horizon problem are achieved through the stress 
energy of a single field, while the primordial perturbations are made up out of the individual 
curvature perturbations of all three fields. 

We will calculate the primordial curvature perturbation after all the particles have de- 
cayed using the 8N formalism and the sudden decay approximation (for a different approach, 
using the ADM formalism see, [IB])- This means that we assume that as soon as the Hubble 
parameter H, becomes equal to the decay rate I\ of a particle, the particle decays instantly 
and transfers all of its energy density (and curvature perturbation) into radiation. In our 
study we will follow the formalism of [TJ5] and |16j : we will obtain a chain of four equations to 
describe the evolution of our model starting from a time when just one particle has decayed 
through to the last particle's decay. Our goal is to find the first and second order components 
of the curvature perturbation at the last decay, expressed in terms of the initial adiabatic 
and isocurvature perturbations: C^,, S a and S x . However, we generalise this method to in- 
clude the linear terms coming from the inflaton, and up to second order terms from the two 
curvatons. Also, different orders of decay for the three particles will be considered, as well 
as situations where one of the curvatons is frozen until after the other two particles have 
decayed. Once we have the expression for £, we can manipulate it in the general form of ^ 
to find the coefficients we need to evaluate f^L- 



3.1 Case 1: t a < t Y < t^, < t a , < t Y , 

u osc Aosc decay u decay Adecay 

In this case, the two curvatons begin oscillating some time after the end of inflation, but 
decay after the inflaton has decayed. We are assuming that the energy density of each particle 
becomes radiation with the particle's decay. When a particle is oscillating, its equation of 
state is the one for pressure- less matter, p = 0, while for a radiation component it is p — |p. 
So, after the inflaton decay but before the first curvaton decay (in other words, while the 
two curvatons are oscillating) the curvature perturbation and energy density of each species 
is 

= (a + -M^) =► P* a = P^~ Ca) (11) 
4 P<i> a 

Co = C. + ^ln(^)^p CTQ =p CTQ e 3 ^) (12) 

Cx = C a + \H^)^ Pxa =p Xa ^-^. (13) 

The index a symbolises the time at which the first curvaton decays and we will use (3 to 
denote the time at which the second curvaton decays, which will be the last decay of the 
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three fields in all the scenarios analysed in this paper. In the equations above, ( a is the total 
curvature perturbation at the first curvaton decay and pi is the average energy density of 
each field i. We can assume that at this point the total curvature perturbation is equal to 



the adiabatic perturbation and so equations ( 12 ) and ( 13 ) become 

5 X 



^ S V -3(C\- U- 



(14) 
(15) 



We don't make use of the above definition for the two isocurvature perturbations immedi- 
ately. We will use equations ( 14 ) and ( [To"]) once we obtain the curvature perturbation at the 
last decay to rewrite it into the form of (l3j). 

Since we are assuming non-linear evolution for the curvatons, p a = \m 
definition c = c(a*), and the parameter d is treated in an analogous way for the second 
curvaton \- On a uniform density hypersurface just before the first curvaton decays, p$ a + 
Pcr a + Px« = Ptotai a • We can divide by p t otai Q and use equations pi]) to (Jl3f to obtain 



2 c 2 , where by 



O ,= 4 (C#.-Ca) _i_ n p 3(C<T-Ca) I n p 3(Cx-C«) _ i 



X" 1 



(16) 



where fL 



Ptotal 



Here we replaced the index <p with 70 in the coefficient Q, since all of 



the inflaton's energy density is now the background radiation. We expand this equation to 
second order which yields three equations. The equation for zeroth order is 



n 



70 a 



A" 



(17) 



which is just an identity; the sum of all the energy densities is the total energy density. The 
first order terms give 



Ca(l 



4Q 



70aS<P(i) 



1) + + 3£l xa ( Xm 



(1) 



where f 10a + + / x 

Ca(2) 



1. The second order part of equation (16) is 

\2 



: ^/^(C^^) — Ca(l)) + 3/<Ta,(Ca-(l) — Ca(l)) + 3/ Xa (Cx(l) _ Ca(l)) Z 
+ Ao a C</>(2) + f<r a C(T{2) + fxcCx(2)' 

After the decay of the first curvaton, a similar analysis gives the equation 



n 



710 



3 4(C 7 i-C Q ) 1 O p3(C x -CcO 



(19) 



(20) 



where now the energy densities of both </> and a together constitute the radiation fluid, 71. 
Expanding this equation to second order gives the following for the zeroth order terms: 



n 



71 a 



(21) 



which we use in the equation that we obtain from the first order part along with (18) to 
finally get the expression below for the incoming radiation perturbation at the decay of x- 

Cn(i) = ^it 1 _ U a - /xdO(i) + Rif<r a (cr(i) + [1 - RiO- ~ /xJ]Cx(i)> (22) 

6 



4-Ov 



where i?i 
eventually get: 



-. We then look at the second order terms and, using equation (19), we 



Cn(2) = — ^(Cti(1) — Ca(l)) 2 + 3(1 — ^l)(Cx(l) ~ Ca(l)) 2 + ^1/700(^(1) ~ Ca(l)Y 
+3-Rl/ CTQ (C(l) — Ca(l)) 2 + 3-Rl/ Xa (Cx(l) _ Ca(l)) 2 + ^i/to«C0(2) 
+-Rl/r Q C<7(2) + [1 - - /x«)]Cx(2)- 



(23) 



Now, considering the decay of the second curvaton x, symbolised by the index (3, we 
obtain 

Onpe^-M + n^e 3 ^-^ = 1, (24) 
valid just before x decays and 

Oyzpe^-M = 1 (25) 

valid right after the last decay. In the equation above, 72 is the radiation fluid after the last 
decay. 

Since all the particles have now decayed into radiation, = C^n) + C8(2) * s the primordial 
curvature perturbation that we wish to evaluate in order to get /jvl- We want to find C/3 in 
terms of the initial, known curvature perturbations of the individual fields. The coefficients 
are going to be functions of the particles' energy densities and of their ratios before the 
second decay and before the third decay. 

We follow the same steps as in the study before and after the decay a: we expand (24) 
and (25) to second order and with appropriate substitutions we find the first and second 
order part of the curvature perturbation at the time of the last decay in terms of C^m; Co-(i) 
and Cxriv The first order component is: 



'(l) - ~ /x/sX 1 - fva - /xJCtf>(i) + #i/<r a (l - /x^)Co-(i) 

+[l-i?i(l-/xJ(l-/x,)]Cx(i) 



(26) 



= C* (1) + \Rif^ ~ /*,)&(!) + \[1 ~ RiO- ~ /xJ(l " fx,)]S x{1) (27) 

where to obtain the last line we used (14) and (15). We can already compare with ^ to 
read off the first three coefficients A, B and T: 



A 
B 

r 



XfiJ 



[1 - R1 (l - f x J(l - f x0 )} 



(28) 
(29) 

(30) 



The equation for the second order part is: 

C/3(2) = 4 ( 1 _ /x^)(C 7 i(l) - C/3(i)) 2 + 3/ X(9 (Cx(i) - C/3(i))' 



(31) 



We won't need any of the results coming from the expansion of (25) but, for completeness, 
we included the calculation in the Appendix. Having obtained equations (26) and (31), we 
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have the expression for the curvature perturbation just after the last decay, which is the one 
we do need: 



C/3 - C/3(i) + C/3( 2 ) 

+ [1 - R^l - / X J(1 - / x ,)]Cx ( i) + 4(1 - / X ,)(C 71(1) - C/^d)) 2 
+3/x^(Cx(i) ~ Cs(i)) 2 + (1 ~ fxp)C r yi(2) + fxp(x(2)- (32) 



To be able to read off the necessary coefficients using (J3J), we just need to write (32) in 
terms of the quantities C<£m; ^(l) an d ^xay ^° acn i eve this we make use of equations (14), 



(15), (22), (23) and (|26j). We can ignore the second order term of the inflaton's curvature 
perturbation C<£(2)' as i n 02]' ^ut we will keep the linear term so = C<t> w - As for the 
curvatons, we consider first and second order terms, d = + |Ci (2) 5 where i = a,x , 
keeping in mind equation §2§. 

Below, we list the remaining 6 coefficients: 

A = ^[-2 + 2R 1 -R 1 f tTa (l-f X0 )-2R 1 f Xa (l-f x2 ) + f xp -2R 1 f X0 

+RiUf^ + ((i - RiO- - /*,)) + f XP )%\ (33) 

+3(1 - R 1 )f 2 aa + mflfy 0a + SR 1 f aa (l - u a y 
3 cc" 

+3RiflJ Xa -- 2 RiUSi-^)] (34) 

Z = ^[^/l/; + 3^/1(1 - - 4(1 - ^) 2 (1 - / J 2 

+3(1 - R 1 )(l - f Xa ) 2 + ARJ^Jlf^ + mfa a & 

+ \^kAl- f x J{l-%)]-\{l-%) (35) 

cc" 

H = -R lfaa (l - f Xf> ) + Rl f aa — (36) 

dd" 

= _ 2 + 2R 1 -2R 1 f Xa (l-f X0 ) + f Xp -2R 1 f x? + ((l-R 1 (l-f Xa )) + f X0 )-- T (37) 



d' 

I = ^^[-8 J R 1 /. Q (l-/xJ-6^ Q (l-/xJ/x,(l-/x^ 



+8(1 - i?!) 2 /^(l - /.J - 6(1 - fli)/„ a (l - /xJ 
+8R 1 f 70a f (Ta f Xa - QRifa a f Xa (^ - fa a ) 

-SRiUJx^ ~ fxJ\ (38) 

We now have everything we need to evaluate the non-linearity parameter for the given 
configuration, making use of equation Q. Apart from the nine coefficients that we have 
evaluated above, we need to use the Planck mass Mpi = 2.4 x lO^GeV^, the first slow-roll 
parameter e* = lMp ; (y) 2 as well as c and d. In the case of a quadratic potential for both 
the curvatons, it is simply c = a and d = \- 
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3.1.1 Second Consistency Check - f NL for a single curvaton. 



Having derived these parameters for this case, we are now in a position to be able to make 
two more consistency checks which compare our results to previous work. We can assume 
that the second curvaton \ is not present and that the contribution to the perturbation from 
the radiation field (i.e. the contribution from the inflaton) is equal to zero £ 7 = 0. We have 
A = r = A = Z = H = 6 = I = 0, and (|9b becomes 



( i+ £)-§-i'- (39) 

which is the same as the single curvaton expression found in [19] (and in |20j for cc"/c' 2 = 0). 

Note that the fractional contributions to the density of the radiation coming from the 
decay (e.g. f aa ) depend in the simple harmonic case upon the mass, the decay rate and the 
expectation value during inflation of the curvaton field. In both of the situations we looked 



at above (here and in section 2.2.1 ), the non-gaussianity does not therefore depend upon the 
expectation value of the curvaton field explicitly, although in the case of two curvatons, it 
does depend upon the ratio of the expectation values. 



3.1.2 Third Consistency Check -fNL for mixed perturbations from one inflaton 
and one curvaton 

If we assume that only the inflaton <p and a single curvaton a are present, we immediately 
see that r = Z = @ = I = and that the rest of the coefficients take much simplified forms. 
In the end, our expression for f^L is reduced to: 

5 32 

fc = 6|9ff ; + mit.fi? x {9g - " <r M "^<- 3 + if - + 2/ °-> (40 > 

+4M 2 pl t,a 2 ,UA2f a A-3 + 2/„„ + PJ - /„„(-3 + 4/„„ + 2fiJ) (41) 

In the situation where there are two curvatons and no inflaton, the non-gaussianity depends 
upon the ratio of the two expectation values of the curvatons. In this case of one inflaton 
and one curvaton, the expression for fNL depends on the curvaton expectation value; we 
study it for a* = 10 14 GeV>* = 10 15 GeV, a* = 10 w GeV and a, = 10 17 GeV in plots gaud 

El 

For the first three cases we notice that the non-gaussianity increases as f aa becomes 
smaller, but not indefinitely. Instead, f^L has a maximum value of fNL ma x = 1298, fNL max = 
128.8 and fNLmax = H-9 respectively. So if we were to observe a non-gaussianity of fNL = 50 
which we were to try to interpret in terms of mixed perturbations from a single inflaton 
and a single curvaton both with harmonic potentials, there would be two combinations of 
mass/decay time which would correspond to the same fNL- For er* = 10 17 GeV we see in 
Figure [2] that the non-gaussianity starts from values around -1 and as f Ua becomes smaller 
it increases and reaches a maximum of fNLmax = 0-83 around which it then stays constant. 

This is the same model as the one studied recently in [17] ; we see the same capacity 
of the model to produce large values of /jvx and the same increase of the non-gaussianity 
parameter when the curvaton expectation value becomes smaller. We also find that Jnl has 
a maximum possible value and it in fact returns to the region allowed by observations after 
f a is smaller than 2 x 10~ 3 for a = 10 15 GeV. Our Figures [l|and @ are in agreement with 
their Fig. 2. 
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0.01 0.02 0.05 0.10 0.20 0.50 1.00 0.10 0.15 0.20 0.30 0.50 0.70 1.00 

Figure 2: The non-gaussianity f NL for f aa G [0,1], for cr* = 10 16 CreV on the left and 
a* = \Q l7 GeV on the right in the 1 inflaton-1 curvaton model. 
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Figure 3: Contours for fNL(fa a , / Xj3 ) as predicted in a simplified version of Case 1 for cr* = 
10 u GeV, f aa G [0, 1] and f Xp G [0, 1] on the left. In the second plot, cr* = 10 w GeV while 
the range is f Ba G [10~ 3 , 1] and f Xfj G [10~ 3 , 1]. 



3.1.3 Case 1 with simultaneous decay 

We would like to put some numbers into our equations to see what they are telling us. In this 
section, we will evaluate the f^L under the simplifying assumption that the second and third 
decay of case 1 happen simultaneously. Furthermore, we assume that both curvatons have a 
perfectly quadratic potential so that cc" /c' 2 = and dd"/d' 2 = 0. To completely simplify the 
situation we assume that the two curvatons have equal expectation values at horizon crossing, 
o* = x*- m our plots for f NL we use the reduced Planck mass M P[ = 2.4 x 10 18 GeV and 
the slow roll parameter e* 0.02. 



After the inflaton (j) decays, the two curvatons decay simultaneously, Q 



This 



means that t a = tp, but to remain consistent with our notation we will express the curva- 
ture perturbation coefficients in terms of / CTq and f x „. In reference |16j . we find that for 
simultaneous decay of the two curvatons the following equation holds: 



fxa 



(42) 



which immediately significantly simplifies the nine coefficients (A,B,r...) in the expression 
for C @. The expression above is a direct result of the fact that we can consider the value 
of f Xa as constrained by the given values of f <Ja and f Xfj . Since / 7oQ + f aa + f Xa = 1 we can 
write the inequality f x < 1 — f a and this eventually leads to /. ^ ri - - 



< ( 1 + L ^ L )f XB which is 



valid for all configurations. We use the simplified coefficients from equations (A-6) to (A-14) 
and equation (|9) to obtain Figure |3]for f aa G [0, 1] and f Xp G [0, 1]. 

For cr* = l(r A GeV, we find that the maximum value for the non-linearity parameter is 
fNL m ax = 1298. This value of Jnl is obtained when f x „ goes to zero (or alternatively f aa —> 
since this case is symmetric), while the contribution of the other curvaton is ~ 5 x 10 -4 . 
This is precisely the same maximum value of /a/x as the one in Figure [T] where we studied 
a model with one inflaton and one curvaton, as would be expected. It's the same in the case 



of o* = 10 GeV ; we find that the maximum value of the non-gaussianity is / 



NL r 



11.9 
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and is identical to the result we got in the 1 inflaton 1 curvaton case. This maximum value 
is not observed in the non-gaussianity produced in the two curvaton scenario as studied in 
[16]. The only difference between our Case 1 and the model in [16J is the inclusion of the 
inffaton's contribution to the perturbations in our work, which creates an upper boundary 
for f NL . 

It is evident in all our plots that there are many different combinations of f aa and f x „ 
that can give the same prediction for non-gaussianity. On the other hand, we want to 
emphasise with Figure [4j that even when we keep the total contribution of the curvatons 
to the perturbation constant (red dashed lines), the /nl can still change significantly. 




0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 0.0012 
ft. 



Figure 4: The area of large f^L in case 1 for a* = 10 GeV. The red dashed lines correspond 
to fa a + f X B = constant. 



3.2 Case 2: t a < t Y < t a , < t<k, < t Y , 

u osc Aosc u decay H^decay \d~ 



ecay 



For the second case, we assume that one of the curvatons decays early, before the inflaton. 
Following the same procedure as in the first case we write down an equation valid just before 
the inflaton decays (decay 1): 

fi 7 oae 4(CCT - Ca) + tt^e 3 ^-^ + Q^e 3 ^-^ = 1. (43) 

Here, the energy density of the curvaton a, which has decayed, forms the background ra- 
diation 7q, while the other two particles are still oscillating so their energy density evolves 
like matter. This means that in what follows in this section, Q loa = Vt aa and f l0a = f aa = 
1 — f<i> a ~ /fe- Once the inflaton decays as well, we have 

7ia e 4(( ^ (J + Q^e 3 ^-^ = 1 (44) 

where the energy densities of <fi and o have been converted to radiation, and we see that the 
system is now the same as in the first case: the inflaton and one curvaton have decayed into 
radiation, while the second curvaton oscillates. This means that in this case the equations 
valid before and after the last decay are simply (24) and (25) respectively. To obtain the 
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coefficients present in the non-gaussianity expression ([9]), we need the first and second order 
components of the curvature perturbation at the last decay. Just like we did in the previous 
case we expand the equations above to second order, to eventually get 



+ [1 - R^l - / X J(1 - f x ,)K x{1) + 4(1 - / v )(Cn ( i) - (p m Y 



+3/x/3(Cx(l) Csm) + (1 /x,s)C7i(2) + /x/sC 



>X(2)> 



(45) 



which we can rewrite solely in terms of C<Pimi ^o-(i) an d ^xm ^° reac ^ °^ the curvature 



perturbation coefficients. The list can be found in equations (A-15)-(A-23) in the Appendix. 



3.2.1 Case 2 with simultaneous decay 

Again, we would like to evaluate our equations to find out what new information they are 
giving us. 

In this case we assume that after curvaton a decays, the inflaton and the other curvaton 
decay at the same time. The coefficients in the curvature perturbation expression written 
in terms of f aa and f x can be obtained by placin g the r elat ion / 7 (to = f aa = 1 - f$ a - f Xa 



into equations (A-15)-(A-23) to produce equations (A-24) to (A-32). We then create contour 
plots of our expression for f^i in Figure [5j The maximum value for f^L is /vL max = 2162.5 
for er* = 10 1A GeV. This is significantly larger than the maximum f^i value of Case 1, which 
shows that the order of the particles' decay can have a big effect on non-gaussianity. 
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Figure 5: f^L in Case 2 when the second and third decay happen simultaneously, in terms of 
f Ua and f x „. The plot on the left is done for cr* = lQ 1A GeV with f aa e [0, 1] and f x „ G [0, 1], 
while the plot on the right for a* = 10 16 GeV with f aa e [1(T 3 , 1] and f x e [HT 3 , 1]. 



We see that the prediction for f^L is significantly smaller for larger cr*. In the case where 



cr* = 10 GeV we find the maximum value / 



NL r 



19.2. 



13 



3.3 Cj3.sc 3i t(j <C t-y <C tn- . <C t^j ^nij 

u osc Aosc u decay Adecay *r decay 

In the third scenario, we consider a late decaying inflaton, decaying after both the curvatons 
have oscillated and decayed. The decay of the curvaton a creates the radiation fluid present 
just before the decay of x- 

O 70Q e 4(c '- CD) + Q xa e 3 ^-^ + Q^e 3 ^'^ = 1. (46) 

After the second curvaton decays as well, we have: 

fi 71Q e 4(C71 ~ Ca) + £V e3(C *~ Ca) = L ( 47 ) 

We can also write the following two equations valid before and after the last decay: 

Sl^pe^-n-M+n^e^-tp) = 1 (48) 
Q^e^-^ = 1. (49) 

These are the four main equations for this case and expanding them to second order even- 
tually results in the expression for the primordial curvature perturbation at the time of last 
decay 

C/3 = [1 - Rli 1 ~ /</>J(l - /^)]C^(l) + -Rl/<T„(1 - /^)Ccr(l) 

+-Ri/ Xq ( 1 _ /^)Cx(i) + 4 /^(C 7 i(i) - C/?(i)) 2 

+3(1 - /^)(C<^(1) - C/3(i)) 2 + /^C 7 l(2) + (! ~ f<t>p)(t(2) ( 50 ) 

We rewrite this to match the form of equation ([3]). We list the resulting nine coefficients in 
the Appendix. 

Case 3, in the simultaneous decay limit is equivalent to Case 2 for b -h- <p. 



4 Discussion 

In this work, we have studied the non-gaussianity predicted in models of inflation with 1 
inflaton and 2 curvatons. We derived a general expression for Jnl from the bi-spectrum of a 
general curvature perturbation (. We saw that our general expression simplifies to reproduce 
other expressions for f^L i n the literature under the relevant simplifications. 

We then analytically derived the necessary coefficients for specific scenarios in which we 
varied the sequence of the decays of the three particle species. For our main three cases, we 
then simplified our fisri expression under the assumption that the last two decays happen at 
the same time. 

For all of the calculations we performed in this paper, we assumed a harmonic potential 
for the inflaton of the form m 2 2 and we assumed therefore a slow roll parameter 60 efolds 
before the end of inflation of e ~ 0.02. We made no assumption for the potential of the 
curvatons in our general calculations; for the simplified cases we used the quadratic form. 

In the simultaneous decay limit we studied the features of the predicted /jvi, and found 
that the produced non-gaussianity is bounded. When the contribution of the curvaton to 
the final perturbation is of order unity, /atl is slightly negative but becomes more positive 
as f 0a and f X0 become smaller. The maximum value of /ml is set in this situation by the 
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expectation value <r* of the curvaton during inflation and is smaller when a* is larger. For 
expectation values smaller than around 10 16 GeV f^L grows beyond the current observational 
constraints if the relative contribution of the curvaton has the right value. If the expectation 
value of the curvaton is 10 17 GeV then we find that the maximum non-gaussianity is around 
Jnl ~ 0.8 which is not in the range that will be detectable by Planck. 

We find that /tvl can take very large values, depending upon the expectation value of the 
curvaton during inflation but that it does not grow indefinitely. After it reaches a maximum 
it decreases and returns to within the range allowed by current observations. 

A particular value of /nl therefore corresponds to two different sets of parameters, even 
in the situation where there is a single inflaton and a single curvaton and more information 
is required. 

The predicted /nl in all of the three field models in this paper is quite a complicated 
function of the relative contributions of the two curvatons. If, for example, we take the 
situation where 80% of the total curvature perturbation comes from the inflaton and 20% 
comes from the combined effect of the two curvatons, there are situations where having all of 
the curvaton contribution to the pertubation coming from one of the two curvaton fields will 
maximise Jnl, arid there are situations where having 10% of the perturbation coming from 
each curvaton will maximise Jnl, depending upon the expectation values of the curvaton 
fields. There are therefore in general several combinations of field parameters that give the 
same non-gaussianity and, alternatively, even when these parameters change in a way that 
leaves the total curvaton contribution to ( constant there can be big changes in the value of 
J'nl- 

Of course it would be interesting therefore to investigate the trispectrum to see if the 
parameter g^L could contribute to distinguishing between different models. It is also quite 
clear that the machinery used in this paper leads to rather a lot of long equations. It would 
be nice to come up with a way of simplifying the calculations such that one could make more 
general statements about N curvatons. 
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A Appendix 
A.l Case 1 



Expanding equation (25) up to second order yields: 



1 - ^72/3 [1 + 4 (C72(l) + oC 7 2( 2 ) ~ Cflm - oCS( 



16 

+y(C72 ( i)-C/3 (2) ) 2 ]- (A-l) 
Zeroth order is just 

^72/3 = 1 (A-2) 

which is an identity, as after the decay of the second curvaton all three particles have decayed 
into radiation. The first order terms give simply 

0(1) = Ct2(i); (A-3) 
while the second order terms give the equation 

^72(2) = -2 (C72(1) - C/3( 2 )) 2 + C/3( 2 )- ( A " 4 ) 

Combining these two equations we get 

£72(2) = Cfl(2)i (A-5) 
which is consistent since all that's left after all the particles have decayed is radiation. 
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A. 1.1 Case 1 with simultaneous decay 

Below we list the nine curvature perturbation coefficients for Case 1 when we assume that 
the two curvatons decay at the same time. The coefficients are given in terms of j Ga and 

fxp- 



A = 
B = 


1 

f<J a 

3 


(A-6) 
(A-7) 


r = 


g (3 + fa a )fxD 


(A-8) 


A = 


^[-3/^-/^(3 + 2/^,)] 

Z 


(A-9) 


E = 


-^4 [-9 + 6/ CTa (2 + f xp ) + 2/* a (3 + f X0 )] 


(A-10) 


Z = 
H = 
= 


-^/x,[-81 + 6(2 + / CT J(3 + f aa ) 2 f xp + 2(3 + f aa ff xp ] 
-^(3 + 2/,J/ X/3 


(A-11) 
(A-12) 
(A-13) 


I = 


-^/. Q (3 + / CT J/ X ,[3(2 + f X0 ) + /<7 a (3 + / X/3 )] 


(A-14) 
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.2 Case 2 

ie nine coefficients needed in the /nl expresion ^ are listed below: 

A = 1 (A-15) 

B = 1^(1-/^) (A-16) 

= |[l-«i(l-/x-)(l-/ w )l ( A - 17 ) 

3 cc" 
= -[-1 + R, - ^(1 - / X/3 ) - it^l - / Xfl ) - R 1 f X0 + R 1 f aa (l - — ) 

+(1 " Ri{l ~ /xJ)(l " ^) + /x,^] (A-18) 

= ^.(1 " / w )[5fli + 2i?i/i + 2/ <r „(-l - 4tf? + 3i??/ x , + 
1 „„// 



r 

A 
E 



Z = 



+l R if°^ (A-19) 
9 



9 {-% + ^(l-W[-5 + 13i?i-8^ 

+6^/ x , + + / Xa (4 - 17i2i + 16iJ? - l2R\f Xp - 4R?/J,)] 

+2&(-l + ZR X - AR{ + R 1 f aa + 3R 2 J xp + Rff^))} 

+ (A-20) 
cc" 

H = - Rl f aa (i - f^) + R^— (A-21) 

6 = + + + (A-22) 

I = ^(1-/^)11-5^ + 4^-3^ 

-R\f X(i + /x«(-l + Ax - + «i/aa + 3i^/ w + J*?&)] (A-23) 
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A. 2.1 Case 2 with simultaneous decay 

With the assumption that the inflaton and the curvaton x decay simultaneously we can 
simplify the coefficients A, B, T, ... and express in the following way: 



A 
B 

r 

A 
E 

Z 
H 


I 



1 

fal 
3 

3 (3 + / x „) 
3 3/ g „ + 4/ x , 
2 (3 + / w ) 

9 2 

1 (4-/ CT J/ x ,[-9 + (21-14/ CT 



2fa a )fxf3.\ 



18 (3 + / x ,) 2 

(4 _ f<T a )fx/3 

(3+7Q 

2 (4 — /<T a ) foafxp 
__ 



A-24) 
A-25) 

A-26) 
A-27) 

A-28) 

A-29) 
A-30) 
A-31) 

A-32) 



(3 + /x,) 
A. 3 Case 3 

The nine coefficients that appear in the curvature perturbation expression for this case are: 



A 
B 

r 

A 
E 
Z 



#l/<Tl(l - / X/3 ) 



-[-1 - 2i2x/i(l - U,) + - + 2i2x/i(l - /,„)(! - 2/ x J] 



(A-33) 
(A-34) 

(A-35) 

(A-36) 



- U,)[5Ri ~ *RifL + 2/a a (-l + 3i2i - 4i?f + 3i^/ 02 + - 6^/^)1(^37) 



18' 

1 r_3_/^ 

9 1 2 



-(l-^ 2 )[-3 + 2i? 1 / CT 3 a +9i? 1 / Xa 



+ 2/' (-1 - AR\ + 3i& 



H 





+2(-l + 2i2i - 4i2? + 3Rlf 4 
+/ ffa (lli2i + 4(-l + J?! - 4i?? + 3?/^ + i??4)/ XQ + 2^4)]} 
-i2i/, a (l - U„)[l + 2/ CT 2 a - 2/ ffa (l - 2/ x J] 
-l + fli/^a-W + iZi/^a-W 

-§/a a (i - /^M 4 ^ + ^i/^ 2 + (-1 + fli - 4i?? + m\u P + i??4)/ Xa 

+/ ffa (-l - ±R\ + 3i2?/,„ + + R 1 f Xa )} 



2 <-2 



iRifxJ 

(A-38) 
(A-39) 
(A-40) 

(A-41) 
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A. 4 Case 1*: t a < t&. < t a , <t y < t Y , 

u osc Y decay u decay Aosc Adecay 

This is very similar to the first case we studied in the main body of our paper, except the 
curvaton \ now remains frozen until both <fi and a decayed. Following the same process as 
before, we find the chain of four equations describing our model until the last decay. Before 
a decays, the decay of the inflaton has created the background radiation while x has not 
started oscillating yet and so its energy density hasn't changed from p x : 

Once u decays it becomes part of the radiation too: 

n Tia e 4Kl ^ (a) + Vl XQ = 1. (A-43) 
When x begins moving, but before it decays we have: 

O^e 4 ^"^ + {l^e 3 ^-^ = 1 (A-44) 



which is the same as equation (24), as the sytem is now the same as in case 1. Finally, after 
the last decay the valid equation is just (25). 

We expand these equations to second order to finally find the curvature perturbation at 
last decay: 

C/3 — f~/ 0a (l — fx$)^yo(l) + /°-a(l — /x/s)Cr(l) + fxpCx(l) 

+ 4 ( 1 ~ /x^)(Cti(i) _ C/?(i)) 2 + 3/ X(9 (C X (i) ~~ ^(i)) 2 

+(1 ~~ fxpK-nip) + fxf>Gx(2)- (A-45) 

Once we rewrite this in the form of ^ it results in the following nine coefficients: 



A = 1 
1 

3 



B = Ml -fx,) 



r 



fxp 
3 



3 rr" dd" 

a = ^ - ^)[-(i - ufi - UaO- -^)-a-^)i 

1 cc" 
E = [8 + 5f 2 aa + f aa (-13 + 6f X0 + 2fl) -3(1-—)] 

1 1 dd" 

Z = gfxA 3 + fxp) ~ Q^XfA 1 ~ -jp) 

Cc" 2 
H = ~fa a (l - - (1 - fa a )f aa 

9 = 

2 
9 



(A- 


•46) 


(A- 


47) 


(A- 


48) 


(A- 


49) 


(A- 


•50) 


(A 


-51) 


(A 


52) 


(A 


-53) 


(A 


54) 
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A. 5 Case 2*: t a <t a . < U. < t y < t Y . 

u osc u decay *r decay Aosc /idecay 

In the final scenario that we study, the curvaton x still remains frozen until the other 
two particles have decayed. The inflaton is no longer the first one to decay, making this a 
modified version of case 2. The equation below, is valid just before decays, so the radiation 
component has already been created by the a decay: 

Sl^e^-C*) + n aa e 3 «°-^ + Q xa = 1, (A-55) 

and once the inflaton decays as well we have the equation: 

n 7ia e 4 ^- Ca > + n xa = 1. (A-56) 

Eventually, x begins to oscillate and so its energy density will evolve like pressureless matter 
and equation (24) is valid to describe the era until just before it decays. After x decays we 
have equation (25) since all that's left is radiation. The primordial curvature perturbation 
we obtain by second order expansion of these equations is: Once we rewrite this in the form 
of (|3f it results in the following nine coefficients: 

A = 1 (A-57) 

B = \U a {l-fxe) (A-58) 

r = \f Xp (A-59) 
3(l-/„) 



A = 2 JX ' 3J [1 ~ - Vl a + Vl a - f Xp + ~ / X ,)f 2 

Z = ^^h^ + 3(1 +441 + ^(1-/^) (A-62) 

rr" 

H = (l-f X0 )[l-2f aa -2fl + 2fl+f aa ^} (A-63) 

© = -«Wx,)(l-^) (A-64) 

I = (1 ~ Q fxp) [~6/ gn (l - f xp )f xp ~ SUJi P ] (A-65) 



+fx^-fx,)-^} (A-60) 
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